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Abstract. In this work we construct compositions of vector processes of the 
form Stf (c 2 L"(t)), t > Q, v £ (0, |] , £ (0, 1], n 6 N, whose distribution is 
related to space-time fractional ra-dimensional telegraph equations. We present 
within a unifying framework the pde connections of ra-dimensional isotropic 
stable processes S^f whose random time is represented by the inverse L v (t), 
t > Q, of the superposition of independent positively-skewed stable processes, 
3V{t) = + (2X)i H%(t), t > 0, (Hf v , H%, independent stable subordi- 

nators). As special cases for n = 1, V = | and (3 = 1 we examine the telegraph 
process T at Brownian time \B\ (Orsingher and Beghin [21]) and establish the 

equality in distribution B (c 2 £5 (t)) l = T(|B(t)|), t > 0. Furthermore the 
iterated Brownian motion (Allouba and Zheng [2]) and the two-dimensional 
motion at finite velocity with a random time are investigated. For all these 
processes we present their counterparts as Brownian motion at delayed stable- 
distributed time. The last section of the paper is devoted to the interplay 
between time-fractional hyperbolic equations and processes defined on the n- 
dimensional Poincare half-space. 
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1. Introduction and preliminaries 

1.1. Introduction. The study of the interplay between fractional equations and 
stochastic processes has began in the middle of the Eighties with the analysis of 
simple time- fractional diffusion equations (see Fujita [11] for a rigorous work on 
this field, or more recently Allouba and Nane [1] , where the compositions of Brow- 
nian sheets with Brownian motions are considered). In some papers the connection 
between fractional diffusion equations and stable processes is explored (see for ex- 
ample Orsingher and Beghin [24]; Zolotarev [29]). The iterated Brownian motion 
has distribution satisfying the following fractional equation 

<92 1 d 2 

— ru(x,t) = zu(x,t), x&R,t>0, (1.1) 

dt? 2s ox 2 

(see for example Allouba and Zheng [2]) and also the fourth-order equation 

= ~u(x,t) + ^^(*), * €R,t > 0, (1.2) 

see DeBlassie [8] (also for an interpretation of the iterated Brownian motion to 
model the motion of a gas in a crack). Zaslavsky [28] has studied the fractional 
kinetic equation (derivatives are meant in the sense of Riemann-Liouville) 

■^p9{x,t) = Lg(x,t)+po(x) _ 0</3<l,zeR, (1.3) 

where po € C°° (M 1 ) is the initial condition and 

Lf = - a2 Tx +Dq d^xr + Dp dx^- (L4) 

For p = q = 1/2, the differential operator (1.4) is symmetric and Saichev and Za- 
slavsky [25] have given the solution to (1.3) in the form g(x, t) = J po(x — y)h(y, t)dy 
where 

/om) = p(x,oh p (-^j r* _1 de, (i.5) 

where p(x,£) is the fundamental soution to 

I - » w 

and the function hg appearing in (1.5) is the law of a positively skewed stable r.v. 
with Laplace transform 

/>oo 

e- st h {t)dt = e~ s \ (1.7) 



Clearly lp(£, t) — jhp [ ly-^r j J £, s 1 is the density of the inverse of H 13 since 

Pr{ff"(f)>£} = Pr{lA0<t}. (1-8) 

Therefore the use of the inverse of subordinators in the solution of fractional equa- 
tions with one time- fractional derivative can be traced back in the papers mentioned 
above and in Baeumer and Meerschaert [3]. 
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When the fractional equation has a telegraph structure, with more than one 
time-fractional derivative involved, that is for v € (0, 1] 

(J^+2\^\u{x,t) = c 2 -^u{x,t), xeR,t>0,X>0,c>0, (1.9) 

the relationship of its solution with the time-changed telegraph processes is exam- 
ined and established in Orsingher and Beghin [21] . The space-fractional telegraph 
equation (with M. Riesz space derivatives) has been considered in Orsingher and 
Zhao [22], while the connection between space- fractional equations and asymmetric 
stable processes has been established in Feller [9] . 

Fractional telegraph equations from the analytic point of view have been studied 
by many authors (see Saxena, Mathai and Haubold [26] for equations with n time 
derivatives). For their solutions have been worked out also numerical techniques 
(see, for example, Momani [20]). Telegraph equations have an extraordinary im- 
portance in electrodynamics (the scalar Maxwell equations are of this type) , in the 
theory of damped vibrations and in probability because they are connected with 
finite velocity random motions. 

In this paper we consider various types of processes obtained by composing 
symmetric stable processes S^f(t), t > 0, < < 1, with the inverse of the 
sum of two independent stable subordinators (instead of one as in Baeumer and 
Meerschaert [3]) say L"(t), t > 0, < v < \. These time-changed processes, 
W„(t) = S 2 f (c 2 ^^)), t > 0, have distributions, w^(x,t), x e R n , t > 0, which 
satisfy telegraph-type space-time fractional equations of the form 

+ 2A^) w£ (x, t) = -c 2 (-Af wl (x, t), xe M", t> 0, c > 0, A > 0, 

(1.10) 

where 0</3<l, < u < \, subject to the initial condition 

wl (x,0) = 6{x). (1.11) 

The fractional Laplacian (— A)^, appearing in (1.10), is defined and analyzed 
in Section 3 below. The fractional derivatives appearing in (1.10) are meant 
in the Dzerbayshan-Caputo sense, that is, for an absolutely continuous function 
/ € L 1 (R) (for fractional calculus consult Kilbas, Srivastava and Trujillo [17]), 

CQV 1 ft dT_rt\ 

-pi — f{t) = tt7 r/ d m-Kv <m,m€N. (1.12) 

dt u w r (m - v) J Q (t - s ) v+1 - m v ' 

Equation (1.10) includes as particular cases all fractional equations studied so far 
(including diffusion equations) and also the main equations of mathematical physics 
as limit cases. Thus the distribution of the composed process S^f (X"(i)), t > 0, 
represents the fundamental solution of the most general n-dimensional time-space 
fractional telegraph equation. We give the general Fourier transform of the solution 
to (1-10) with initial condition (1.11) as 



E u ,i (nt v ) + ( 1 - . A I E uA (r 2 n 



A 2 -c 2 ||C|| 2 V V J\ 2 c 2 1|£" 2/3 



(1.13) 
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where 

n = -A + ^A2- C 2||£|| 2 ' 3 , r 2 = -A - y/x^*Mf*. (1.14) 

and 

fl^oo = E rJx /v ^ >0 ' (1 - 15) 

is the two-parameters Mittag-Leffler function (see, for example, Haubold, Mathai 
and Saxena [15] for a general overview on the Mittag-Leffler functions). Our result 
therefore includes all previous results in a unique framework and sheds an additional 
insight into the literature in this field. 

An important role in our analysis is played by the time change based on the 
process L u (t), t > 0. We consider first the sum of two independent positively 
skewed stable r.v.'s and H%(t), t > 0, < v < \, 

H v {t) = Hl v {t) + (2\)"H%(t), t>0, (1.16) 

whose distribution R v (x,t) is governed by the space fractional equation 

d ( d 2v d v \ 1 

m 6 v (x,t) = -\^ ; +2X—j6 v (x,t), x>0,t>0,0<v<-. (1.17) 

In (1.17) the fractional derivatives must be meant in the Riemann-Liouville sense 
which, for a function /fL 1 (M), is defined as 

d v 1 d m f x f(s) 

— — fix) = — - - — / Htt d s, m - 1 < v < m, m £ N. 

dx v r (m - v) dx m J (x - s ) u+1 ~ m 

(1.18) 

We then take the inverse L v (t), t > 0, to the process H v [t\ t > 0, defined as 

L v (t) = inf {s > : Hf"(s) + (2\)» H%(s) > t} , t > 0, (1.19) 

whose distribution is related to that of H v (t), t > 0, by means of the formula 

Pr {L v {t) < x} = Pr {9T(x) > t} . (1.20) 

The distribution 4,(x, t) of L"(t), t > 0, satisfies the time-fractional telegraph equa- 
tion 

d 2v d v \ d 1 

2\—)C v (x,t)^- — [ v {x,t), x>0,t>0,0<^< -, (1.21) 



where the fractional derivatives appearing in (1.21) are again in the Riemann- 
Liouville sense. We are able to give explicit forms of the Laplace transforms of 
h u (x : t) and C v (x,t) in terms of Mittag-Leffler functions for all values of < v < |. 
For example, for the distribution lv(x,t) of L v (t) we have that, for 7 < A 2 , 

e~~< x C v {x,t)dx = 

1+ , A 1 E v>1 {r x t u ) + I 1 - " ]E Vtl (r a Q , (1.22) 



7, 




where 

ri = -A+ VA 2 -7, r 2 = -A - v^ 2 - 7- (1-23) 
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The distribution Q,{x, t) of L v {t), t > 0, has the general form 
Cu{x,t) = [ hu{x,s) h„(t- s,2\x)ds + 2\ [ l„(2\x, s) h 2u (t - s, x) ds, (1.24) 



where the distributions of H 2u , H v , and that of their inverse processes l? v and 
L v appear. For our analysis it is relevant to obtain the distributions of J{^(t), 
t > 0, and L^{t), t > 0. We also obtain explicitely the distributions of Hz{t) and 
Hi(t), t > 0, and also of their inverses £s(i) and Lz(t), t > 0, in terms of Airy 
functions. By means of the convolutions of these distributions we arrive at the 
following cumbersome density of the random time t > 0, 




(1.25) 



For n = 1, /? = 1 and v = 1 in (1.10), we get the telegraph equation which is 
satisfied by the distribution of the one-dimensional telegraph process 

T(t) =V(0) f (-l) N ^ds, t>0, (1.26) 



where N(t), t > is an homogeneous Poisson process, with parameter A > 0, 
independent from the symmetric r.v. V(0) (with values ±c). Properties of this 
process (including first-passage time distributions) are studied in Foong and Kanno 
[10] and a telegraph process with random velocities has been recently considered 
by Stadje and Zacks [27]. 

For n = 1, /3 = 1 and v = i the special equation 

(l+ 2A ^l) w |(^*) = c*£,w\(x,t), xeR,t>0, 
w\(x,0) = 5(x), 

v 2 

has solution coinciding with the distribution of T(\B(t)\), t > 0, where |£>(t)|, 
t > 0, is a reflecting Brownian motion independent from T (see Orsingher and 
Beghin [21]). For A — > oo, c — > oo, in such a way that — > 1 the fractional 
diffusion equation (1.1) is obtained from (1.27) and the composition T {\B(t)\), t > 
0, converges in distribution to the iterated Brownian motion. Our result, specialized 
to this particular case gives the following unexpected equality in distribution 

T(\B(t)\) 'l= w B^dit)), t>0, (1.28) 

where 

Pr {B (c 2 L»(t)) € dx\ = — f - 1 e-ik~^ ( ^— + da, 



and 



/>oo 

Pv{T(\B(t)\) e dx} = / Pr {T(s) G dx} Pr{|S(t) | G ds}. 
Jo 



(1.29) 
(1.30) 
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The absolutely continuous component of the distribution of the telegraph process 
T(t), t > 0, reads 



where \x\ < ct, t > 0, c > 0, and 



(1.31) 



= £ (If w- c 1 - 32 ) 



fc=0 



(fc!) 



For n = 2, /3 = 1 and v = 1, equation (1.10) coincides with that of damped 
planar vibrations (we call it planar telegraph equation) and governs the vertical 
oscillations of thin deformable structures. The solution to 



' ( |£ + 2A|) r(x, y, t) = c 2 ( + ^) r(s, y, i), z 2 + y 2 < c 2 t 2 , t > 0, 
r(a;,2/,0) = 5(x,y), 
r t (x,y,Q) = 0, 

(1.33) 

corresponds to the distribution r(x, y, t) of the vector T(t) = (X(t), Y(t)) related to 
a planar motion described in Orsingher and De Gregorio [23] . This random motion 
T(t), t > 0, is performed at finite velocity c, possesses sample paths composed by 
segments whose orientation is uniform in (0, 27r), and with changes of direction at 
Poisson times. The distribution r(x,y,t) of T(t), t > 0, is concentrated inside a 
circle C ct of radius ct and has an absolutely continuous component which reads 



r(x, y, t) = - (z, y) g C rf , t > 0. (1.34) 

2?rc yJcH 2 - (x 2 + y 2 ) 

If no Poisson event occurs, the moving particle reaches the boundary dC c t of C c t 
with probability e~ A *. The vector process T(t), t > 0, taken at a random time 
represented by a reflecting Brownian motion, | J3 (i) | , has distribution 

/>oo 

q(x,y,t) = / Pr {X{t) eds,Y(t) eds}Pr{\B{t)\eds} (1.35) 
Jo 

which satisfies the fractional equation 

(0 d" \ ( Q 2 d 2 \ 

at + 2A atI ) q{x > y ' t) = <? \d^ + d^) q{x > y ' <} ' v) e R2 ^ > °- 

(1.36) 

However, the distribution of £?2 (c 2 £5(t)^ ; t > 0, does not coincide with (1.35) 

(£?2 is a two dimensional Brownian motion). In this case the role of T(t), t > 0, in 
(1.28) is here played by a process which is a slight modification of T(i), t > 0. We 
take the planar process with law 



. \e~ xt 
r(x,y,t) = 



^c 2 i 2 - (a; 2 + y 2 ) 



:r 2 +y 2 < c 2 t 2 ,t > 0, 
(1.37) 
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which also solves equation (1.33). The process with distribution 

r°° r i 95 

q(x,y,t) = / t(x,y,s) Pr {\B(t)\ e ds} + — Pr e rfs} 

JO [ ■ZAOtz 

t(x, y, s) + JU(ar, y, a)j Pv{\B(t)\ e ds} , (1.38) 

has the same law ol a planar Brownian motion at the time L^(t), t > 0. The 
process T(t), t > 0, possessing distribution (1.37) is obtained from T(t), t > 0, by 
disregarding displacements started off by even-order Poisson events. 

The last section of the paper is concerned with random motions on the hyper- 
bolic Poincare half-space, H™ = {x,y : x G > 0}, whose distributions are 
governed by fractional equations of the form 

PrMt) = % (sinh"" 1 ^ , r, > 0,t > 

(1.39) 



a 2 " _j_ 9 \ 



IX (?7,0) = 5(J7), 

for < ^ < | and n G N. The corresponding kernel 

<(*?,*) = . , Li *7>0,i>0, (1.40) 

sinn t] 

solves instead the fractional equations 
f(j^+2Af;)<M) = iis^^^h"- 1 ^^,*)), 77>0,i>0, 

l<M) = 5(7,). 

(1.41) 

The process t > 0, in HP which possesses distribution p£(x,t) solving (1.39) 

is obtained by means of the composition 



^{t)=B^{L"{t)) ) *>0, (1.42) 

where -B^ p is the hyperbolic Brownian motion in EP. The hyperbolic Brownian 
motion has been introduced in the plane by Gertsenshtein and Vasiliev [12] and in 
H 3 by Karpelevich, Tutubalin and Shur [16], in 1959. In successive papers many 
properties of the hyperbolic Brownian motions have been explored (see for example 
Getoor [13]; Gruet [14]; Lao and Orsingher [18]; Matsumoto and Yor [19]). The 
relationship between kernels in H 2 and H 3 and kernels in higher-order spaces is 
represented by Millson formula 

e -nt Q 

k„+2{v,t) = - „ r-r ttMM), 7/ > 0, i > 0, n G N. (1.43) 
2tt smn 77 or] 

Since p^ p and k$ are considerably simpler than p^ p and fc 2 we give explicit expres- 
sions for the distribution 

1 , Ansinhji f* e~ s a 2 , 2 >, 2 / s \ , . „ 

ft 1 (»?,*) = — 4 / 3 e~— *• +2 ds, (1.44) 

where 77 > and i > 0. This distribution solves the fractional- hyperbolic telegraph 
equation (1.39), for v ~ \ and n = 3. 
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1.2. Notations. For the reader convenience we list below the main notations used 
throughout the paper. 

• Sf{t) = (sf{t),S 2 2 p {t),--- ,Sf{t)), t > 0, < < 1, n e N is a 

isotropic stable n-dimensional process with law vp (x,t), x G M. n , t > 0. 

• H v (t), t > 0, < v < 1, is a totally positively-skewed stable process (stable 
subordinator), with law h v (x,t), x > 0, t > 0. 

• L v {t), t > 0, is the inverse of H v {t), t > 0, and has law l v (x,t), x > 0, 
f > 0. 

• 9C{t) = Hl y (t) + (2A) " H%(t), t > 0, is the sum of two independent stable 
subordinators and has law ^{x^t), x > 0, t > 0. 

• L v (t),t > 0, is the inverse of H v {i), t > and possesses distribution C u (x, t), 
x>0,t>0. 

• T(t), t > 0, is a telegraph process with parameters c > and A > and 
\&w px(x,t), —ct < x < ct, t > 0. 

• W n (t) = Sf (c 2 L u (t)), t > 0, has law (x,t), xeR n ,t>0. 

• <W{t) =T(\B(t)\), t > 0, has distribution w(x,t), x G E, i > 0. 

• T(t), i > 0, is the planar process with infinite directions, parameters c, A > 
and law r(x,y,t), (x,y) G C ct = {{x,y) e M 2 : x 2 + y 2 < c 2 t 2 }, t > 0. 

• X(i), t > 0, is the planar process with infinite directions, parameters c, A > 
and law t(x,y,t), (x,y) G C ct = {0,y) el 2 :i 2 +j 2 < c 2 i 2 }, £ > 0, 
constructed by disregading displacements started off only by even-labelled 
Poisson events. 

• Q(t) =T(\B(t)\), t > 0, has law q(x,y,t), (x,y) G M 2 , t > 0. 

• S„ p (t), t > 0, is the n-dimensional hyperbolic Brownian motion in H 71 = 
{{x,y) : x G M n_1 ,y > 0} and has law p^ p (?7,i), ?y > 0, t > with kernel 
k„(r),t), i] > 0, t > 0. 

• T^(t) = B^(L u (t)), t > 0, has distribution f^foi), 77 > 0, i > and 
kernel K^(rj,t), tj > 0, t > 0. 

• By / we denote the Laplace transform of the function / and by / we denote 
its Fourier transform. 

1.3. Preliminaries. Let us consider a stable process ^(i), t > 0, < v < 2, 
v 7^ 1, with characteristic function 

fl5gt£S"(t) _ e -CT|f|"t(l-i9sign(f) tan ^) (1-45) 

where 9 £ [ — 1, 1] is the skewness parameter and 

■kv 

a = cos — . (1.46) 

For 8 — 1 the distribution corresponding to (1-45) is totally positively skewed and 
for 6 = —1 is totally negatively skewed. The stable process with stationary and 
independent increments, totally positively skewed will be denoted as H v (t) : t > 0. 
We note that the density h v (x,t), of H u (t), is zero at x — as the following 
calculation show 

1 r°° 1 r°° 



fc„(0,t) = — / Ee* H "Wdg = — / e-^ISr^i-ite" 1 ^)^ 

/■O 

r|$rt(l-itan i 5 L )^ _|_ / e -<r|f|"t(l+itan 



1 

27 



e 
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1 

2^ 
1 

2^ 
cos ■ 







--1 



e 2 dz 

z\b-i l 



-i l 
t 



2 dz 



t 



dz = 0. 



The positively skewed stable r.v. H u (t) has a:-Laplace transform 
MM) = Ee-^ ff "W = e - *"", < v < 1, 
and therefore Fourier transform 



(1.47) 



(1- 



(1.49) 



This shows once again that the skeweness parameter is 9 = 1. 

The probability law /^(rc, t), of H v (t), t > 0, solves the boundary-initial problem 



(1.50) 



(i + J^JMM) = o, x>o,t>o,o<i/<i, 
MM) = o, 

^(x,0) = 

By taking the a;-Laplace transform of the Riemann-Liouville fractional derivative 
appearing in (1.50) we have that 

d v 

Jo dx u 
d f x h v (z,t) 



C 



dx' 



;h u (x,t) 



(M) 



h u (x,t)dx 



1 

r (1 — v) dx J (x — z) v 
1 

T{l-v) J dx z v 



dz 



dx 



d h v (x — z,t) t h v (0,t) 



MM) 

T(1-p)J 
■ h u {0,t)fi 1 '- 1 



-» x x 1 -"- 1 dx + 



1 



1 



r(i-v)J z 



00 dz 



dx 



— / dxe ^ x —h v (x - z,t) 



dx 



T(l-u) 



: h v (0, t)n 



Therefore 



so that 



u-l 



e > ix h u (x 1 t)dx 



J dz 



dx 



h v {x, t)dx 



M MM) = A*"MM) 



' ^h v (fj,,t)+n v h v (fi,t) = 0, fi>0,t>0, 
MM) = 1, 



h v (fj,,t) 



(1.51) 



(1.52) 



(1.53) 



In other words the density of a positively skewed stable r.v. solves the space- 
fractional problem (1.50). 
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We will also deal with the inverse process of H u {t), t > 0, say L v (t), t > 0, for 
which 

Pr{H u (x) > t} = Pr < x} , x > 0,t > 0. (1.54) 

Such a process has non- negative, non-stationary and non-independent increments. 
Furthemore we recall that the law l v (x,t) of L u (t), can be written as 

l v (x,t) = ^W- v ±- v (-^j , x>0,t>0, (1.55) 

where 

oo u 

= E fc ,w! + 5 v ^eM,a>-l ) 6eC, (1.56) 

k=0 ' ^ ' 

is the Wright function, and has Laplace transform 

~ r°° r°° ~\ / r \ 

l v (x,n) = I e-»%(x,t)dt = I e -' d -W- v , 1 - v (--)dt = fi^e^ . 

(1.57) 

2. Sum of stable subordinators, 9{ v (t) = Hf"(t) + (2\)$H%(t) 

For the construction of the vector process W„(i) = S^f (c 2 X K (t)), t > 0, whose 
distribution is driven by the general space-time fractional telegraph equation (1.10), 
we need the sum H v {t\ t > 0, of two independent positively skewed processes. The 
second step consists in constructing the process L u (t), t > 0, inverse to H 1, '(t), 
t > 0. We now start by considering the following sum 

H v [t) = H\ v {t) + (2X)i H%(t), t > 0,0 < v < \, (2.1) 

with Hf v , H%, independent, positively-skewed, stable random variables, A > 0. 
The distribution of H v (t) can be written as 



Hu{x,t) = / h2 V {y,t)h„{x-y,2\t)dy. (2.2) 
Jo 

Taking the double Laplace transform of (2.2), with respect to t and x, we get 



oo />oo 



M7,M) = / e - "* / e- lx h u {x,t)dxdt = / e - "*-* 7 ""^"dt 
■/o ./o Jo 

1 1 



1 



7^ — r% 7" — ri 



7 2l/ + 2A7" + /x 
where, for < /1 < A 2 , 

\r x = -A - v/A 2 - it, 



1 



^2 - ri 



By means of formula 



[r 2 = -A + y 7 A 2 - fi. 

r v 1 

/ e-^x^E^^x") dx = -— 



(2.3) 



(2.4) 



(2.5) 



V 

where E v v {z) is the Mittag-Leffler function defined in (1.15), we can invert the 
x-Laplace transform in (2.3) obtaining, for fj, < A 2 , 

h v (x, fi) = 
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1 



E VM U-X + VA 2 - /,) i") - ((-A - v 7 ^/' I 

i a 



1 d 

-A - ^A 2 - n dx 



A+ a/A 2 - n dx 



Ms 



(2.6) 



Formula (2.6) gives the explicit form of the i-Laplaee transform of /t„(a:, i) in terms 
of Mittag-Leffler functions. In view of formula 



e v , x (-xn = 



■dx, < v < 1, (2.7) 



we have that 



^ (x, fx) 



2\/A^ 



y v 1 sin7r^ dy 



-A + v/A 2 - Jt ^ V ^ (y 2 " + 1 + 2 J/ 1/ cos 7Ti/) 

-anj/f A+a/A 2 -^) " . j 

V v 1 y sm 7tz/ ay 



A + v/A 2 ^ 9a; y 7r(j/ 2 " + l + 2y"cos7rj/) 

dy sin7w 1 



„ 7r(y^ + l + 2y^cos7r^)2 v /A 2 ~^ 



(a- v/a^T) 



1 d 



r2 — ri <9x 



-E 



-ia*(- r 2 ) " 



>2 



-xU"(— ri) f 



?'l 



where W 1 ' is the Lamperti distribution with density 
Pr{« y G(iu} sin.7Ti/ u"^ 1 



du 



u > 0, 



(2.8) 



(2.9) 



7T 1 + u 2v + 2 COS 7Tl/ 

and represents the law of the ratio of two independent stable r.v.'s of the same 
order v. 

Theorem 2.1. The law ft„{x,t) of the process 9i v (i) = Hf u (t) + {2X)» H%(t) solves 
the fractional problem 

§- t h v {x,t) = - (|^7+2A|^) £,(x,f), x>0,t>0,0<^< |, 
/k(0,f) = 0, (2-10) 
/v(z,0) = <5(z). 

The fractional derivatives appearing in (2.10) are intended in the Riemann-Liouville 
sense. 
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Proof. By considering (1.49), we have that the Fourier transform of fk,(x, t) is writ- 
ten as 

__ g — | 2iy cos Txv{\ — i sign(£) tan nu) — 2 At |£| w cos 4p (l — z sign(£) tan ^ ) 



and thus 



at 



>ign(5)\ ^ 2 \t(\£\ e -'-i- Biei,(0 ) 



(2-11) 



. e -t(l«|e-* si8nK) ) 21 



'-2A(|C| £ 
-2At(|5| e -¥ 



iign(5) 



In view of the relationship 

|f| e -Y»w»tt) = _jf 
we have that formula (2.12) can be rewritten as 



0* 



HO -2A Her 



-t(-if) -2Ai(-if)" 



In (1.51) we have shown that 



C 



dx u 



h v (at, t) 



(/'■) 



3 dx" hv ^- x,t ' dx = M h v (fi,t) 



(2.12) 
(2.13) 

(2.14) 
(2.15) 



and thus for a sufficiently good function / we have the following Fourier transform 



T 



dx u 



fix) 



(0 



e-^-^mdx = (2-16) 



In view of (2.16) we have that the Fourier transform of the right-hand side of the 
equation (2.10), equipped with the boundary conditions, is written as 

-h„{x,t)+2\ — h„{x,t) 



T 



dx 2v 



dx u 



(0 



pp.v poo Q.V 

s _( ^raMM) dx -2A / e-^^MM) *c 



(HO 2 H2AHO 1/ )^(C,0 

+ 2A HO") e^O* 1 ""* 



-2At Ifle-TT 



= - (HO 2 " + 2A (-if)") e -*(-^-2A t (-^)^ 



(2.17) 



which coincides with formula (2.14). This is tantamount to saying that the Fourier 
transform h v (£, t) is the solution to 



(2.18) 
□ 



| &k (0 *) = - (HO + 2A (-if)" J ^ (0*) , H M, * > 0, 

lM£,o) = i, 

and this completes the proof. 
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2.1. The inverse process L l ' (t). Let L v (t), t > 0, be the inverse process of ^(t), 
t > 0, as denned in (1.19) for which 

Pr {L"{t) < x} = Pt{!X v (x) > t}, x,t>0, (2.19) 

and let d,{x,t) be the law of L v (t), t > 0. We have the following result. 

Theorem 2.2. The law Ci,{x,t) of the process L v (t), t > 0, solves the time- 
fractional boundary-initial problem 

(U+2Aj£) Ux,t) = —^(u(x,t), x > 0,t > 0,0 < v < §, 
Ux,0) = 5(x), (2.20) 
C(0,i) = r(*i-2iy) + 2/ ^r(*i-i/)' 
and /ias x- Laplace transform which reads, for < 7 < A 2 , 
1 



7^ ^ l(rin+ r7^r" (r2r 



(2.21) 

where 

ri = -A + VA 2 - 7, r 2 = -A - VA 2 - 7- (2-22) 
The fractional derivatives appearing in (2.20) are intended in the Riemann-Liouville 
sense. 

Proof. We first show that the analytical solution to the problem (2.20) has double 
Laplace transform 4, (7, fi) written as 

~ f- 1 + 2A/J"- 1 

By taking the i-Laplace transform of the equation in (2.20) we have that 

H 2 "Z(x,n) + 2\^Z{x, f i) = -^Z(x,n). (2.24) 

By taking into account the boundary condition and performing the x-Laplace trans- 
form of (2.24) we have that 

+ 2\n v ) Z (7, n) = Z (0, //) - 7 Z (7, ft) • (2-25) 
Now, by considering the boundary condition, we get that 

t -2u t -u 

r + 2A— r 

P(l-2v) T(l-v) 

= ^- 1 +2X^- 1 , (2.26) 

and thus 

~ /i 2 "- 1 + 2A^-! 

^ (7>/x) = ^ + 2A^+ 7 - (2 ' 27) 
Now we show that the double Laplace transform of the law 4, (x, t) coincides with 
(2.23). We first recall that 



4,(0,//) = / dte'^X^t) = / dte~^ 



Hu(fJ,,x) = / dte-"%(t,x) = Ee-^W = Ee - "* "(^Ee-^ 2 ^ 
./o 

= ft^ (//, a:) ^ (m, 2Aa;) = e -V"-*2A M " ^ x > Qj ( 228 ) 
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where we used result (1.48). By considering the construction of the process L v (t), 
t > 0, as the inverse process of tt v (t), t > 0, as stated in (2.19), we get 



COM) 



PT{£ v (t) E dx} d 



dx 



(s, x) ds. 

(2.29) 



In view of (2.29), the double Laplace transform of Q,{x, t) can be obtained observing 
that 

d_ " 
dx 



dx, '-■ j dte-f* 
d 



hu{s, x) ds 



dx ( 



ii 



dx 



dte *** / Hu(s,x)ds 



1 

M Jo 

.2i/-l 



_ d ~ 1 

efee ~< x — Itv (x, /j.) = 



da; e" 7 * 



5 



dx 

_ ^ + 2A M 



^ 2l/ + 2A^" + t~' 



(2.30) 

which coincides with (2.23). Now we pass to the derivation of the x-Laplace trans- 



form of 4/ (a?) t)- We can write 

,2u—l i n\,,v—l 



2A/x" 



^+2A/x" + 7 /i"-ri /i'-'-rs (/x w - n) (m" - r 2 ) 



ju" - n \i v - r 2 



fJ> - n \i v — r 2 



2^ 



where 



n = -a + v/a^- 

Now we need the following results 



7: r 2 



~^E Vtl { rj t v ) dt 



,v—l 



-A- VA 2 - 7 . 
j = 1,2, 



(2.31) 
(2.32) 



,2v-l 



M - n 



(2.33) 



Therefore 
£(7,t) = (nf) + (rat") 



27^ 



(nf ) - (r 2 f )] 



Since 



( z ) = zE u x(z) + 



T(l-u) 



(2.34) 
(2.35) 



we have that 

Z(rf,t) = (nf ) + (r 2 f ) 



IT" 



i -a + 7a^ 

2^A^ / 



27A^- 



[r x f (r x f ) - r a f £„,i {r 2 t v )] 



SPACE-TIME FRACTIONAL TELEGRAPH EQUATIONS 



15 



E v ,i [rit v ) + 1 



Va 2 "^ 



E v ,i (r 2 t v ) 



(2.36) 



which coincides with (2.21). 

Now we check that the Laplace transform (2.36) solves the fractional equation 



-7 1(7.*) 



r 



2A- 



r(i-2i/) r(i-v) 

which is the x-Laplace transform of the equation appearing in (2.20). Since 
q2v _ t _2l/ Ca2w _ 

we therefore need to show that 

(CqIv G d v \ ~ ~ 

+ 2A 1 — 4(7,*) = -74,(7, *)■ 



V Si 2 " dt v 



In light of 



St 

Cq2v 



-E v s{r 3 t») = rj E uA ( rj r), j = 1,2, 



<9t 2 <* 



T(l-u) 



(2.37) 

(2.38) 
(2.39) 

(2.40) 

(2.41) 
(2.42) 



we are able to show that (2.21) solves (2.37). We first check result (2.42) as follows, 
for < 2v < 1 



Cq2u 



dt 2 ' 

OO 

E 



tf 2 



fc=0 

z/fc 



r {vk + 1) <9i 2 < 



fe=l 



s vk - y (t- s y 2v d s 



E 



r (i/fc + 1) r (i - 2i/) y 
i 



fc tt r ("*) r(i-2i/) 

oo ±vk—2v 



s vk - l {i- s y 



E 



oo „fe+l jvk—v 



ds 



= — — 



-j r (i/fc - 2i/ + 1) ^ r (i/jfc - v + 1) 



g (^) fc 



^ r (i/fc - 1/ + 1) r(i-i/) 



r 2 (rj-f) + 



* _i/r j 

f(i^)' 



(2.43) 



Therefore 



■Cq2v Cqv 

5^+2A— )4(7,*) 
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Cg 



,2// 



dt 2v 



E v ,i (nt v ) + l 



eg 



2// 



E v ,! (r 2 t v ) 



2A- 



Cqv 



2 

(■• 

f 



+ 2A 2 



(rx*") + 1 - 



A \ C 8 V 



E v< i (r 2 f) 



A 



ri 



r 2 2 E„,i (rat") + 



f(T^) 



A 



ri 1 



A 



r(i-iz) 

(n^,i(rif))+ ( 1 



r 2 £„,i (r 2 i") 



E v ,i (r 2 t u ) (r 2 + 2A) 



(rif ) (rx + 2A) + r 2 1 - 



7 A + y/A 2 - 7 7 ^A 2 - 7 - A 



2 ^-7 



' 7 



1 



A 



In the last steps we used the fact that 



E v ,i {rit v ) + 1 



A 



E v ,i (r 2 f 



1 



A 



1 - 



r 2 t~ v 



^A^7/ ^ ( 1 -^) V 7A 2 ^/r(i- i ,) 



= o, 



and 



7*1 + 2A = -r 2 , r 2 + 2A = — n, = 7. 



(2.44) 
(2.45) 



(2.46) 
□ 



Remark 2.1. The derivation of result (2.21) suggests an alternative proof for the 
Fourier transform (Theorem 2.2 in Orsingher and Beghin [21]) of the law of the 
time-fractional telegraph process. 



Remark 2.2. From (2.31) we get the time Laplace transform of Cu(x,t), for x > 
0,/i > 0,0 < v < \, as 

Since (see formulas (1.55) and (1-57)) 



l v (x,fj) = 

and (see formula (1.53)) 

h u {n,t) 



- /J. I 



1 



f 



dt = n v ~ L e 



V-l -XfJ." 



(2.47) 
(2.48) 



(2.49) 
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we are able to invert (2.47) and we obtain the explicit distribution of the process 
L v {t), t > 0, which reads 

Vx{L v {t) e dx} 



dx 

t rt 
hu (x, s) h v (t — s, 2Xx) ds + 2X / l^(2Xx, s) h 2ll (t — s, x) ds 

Jo 

t 

i / r \ 

ds 



^7 W_ a „,i- 2 » (-—) K{t - s, 2Xx) , 

f f 1 / 2AaA 

+ 2XJ — W- u ,i-u ( — — J h 2v (t - s, x) ds. 



(2.50) 



The densities /i„ and h 2v can be written down in terms of series expansion of stable 
laws (see pag. 245 of Orsingher and Beghin [24]). 

3. n-DIMENSIONAL STABLE LAWS AND FRACTIONAL LAPLACIAN 

Let 

S 2 /(t) = (s 2 \t),S 2 /(t),--- ,S 2 J(t)), i>0,/?6(0,l], (3.1) 
be the isotropic stable n-dimensional process with joint characteristic function 

= e-*ll«ll". (3.2) 

The density corresponding to the characteristic function v^f (£, i) is given by 

vf (x, t) = v 2 J j ( Xl ,x 2 , ...,x n ,t)= — ^ / e-^e-'^l 2 "^. (3.3) 

The equation governing the distribution vff (x,t) of the vector process S^f(t), 
t > 0, is 

(Jt + vf (x,t) = 0, x e R n ,t> 0, (3.4) 

where the fractional negative Laplacian is related to the classical Laplacian by 
means of the following relationships (Bochncr representation, see for example Bal- 
akrishnan [4]; Bochner [6]) 

^ J o °° dX X^ (A-A) -1 A = (jf cto) A dX 

Sin7r/3 Ar(/3) [°° w'-^e-^-^dw = — - a . f°° w^-'e-^dw 



v 'Jo r(i-/3)y 

= -(-A) /3 . (3.5) 

A definition of the fractional negative Laplacian can be given in the space of the 
Fourier transforms as follows 

-(-Afu(x) = { eS{e i +e 2 + ---+gf u{Z)di, (3.6) 

where 

Dom(-A)^ = j«G££,c(R n ) : jfjS (£)| 2 (l + ||£|| 2 ") ^ < oo} . (3.7) 
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An equivalent alternative definition of the n-dimcnsional fractional Laplacian is 

(3.8) 



v> \\x - y\\ 

where the multiplicative constant c(j3,n) must be evaluated in such a way that 



/ e^ x {~Af u(x)dx = ||£|| 2/3 f e* m u(x)dx. 



(3.9) 



Let us focus our attention on the one-dimensional case of (3.8). In this case we 
have that, for < 2/3 < 1, 

d 2 ^P 



dx 2 



u{x) = c(/3,l)P.V. 



u(x) - u(y) 
\x - y\ 



1+2/3 



= c(/3,l) lim 



= c(/3,l) lim 

_ r(i - 2/3) 



2/3 



0+e 

c(/3, 1) 



e w(a;) — u(x — z) 

u(x) — u(a; + z) 

1 d 
r(l-2/3)cfc 



<iz + 



u(x) — u(x — z) 



0+e 

oo 



u(x) — u(x — z) 



dz 



rl + 2/3 



dz 



u(z) dz 
{x ~ z)W 



u{z) dz 

{z~xyp 



(3.10) 



where in the intermediate steps, we considered the relation between the Marchaud 
and the Weyl fractional derivatives. By setting 

2/3 



c(/3, 1) = 



2 T (1-2/3) cosA3tt' 



(3.11) 



we have that, for < 2/3 < 1, 



dx 2 
1 



2 cos (3ir 
1 



u{x) 



1 



x u(z) dz 



r(l - 2/3) dx (x - z) 2 ? r(l - 2/3) cfo ^ (z - x) 



u(z) dz 



)2ri 



u(z) 



2 cos /3tt T(l - 2/3) cfa; 7.^ |x - z\ 

Q2H 



\2n 



dz 



d 2 ? 

d\x\ 2 ? 



u(x), 



(3.12) 



where represents the Riesz operator. 

Remark 3.1. We notice that, for < 2/3 < 1 

Q2P 



T 



_d\x\W 

This is due to the calculation 

Q2fi 



u[x) 



(0 = -iei 2/3 «(e) 



(3.13) 



J" 



1 



(0 



1 



2cos/37i-r (1-2/3) 
2cos/37i-r(l - 2/3) 



die** 

1 

dx e^ x 



d 
dx 



u(z) dz 



u{z) dz 
(x - zf !i dx J x (z-x) 



x u(z) dz 

^8 



(a; — z) Jie (z — x) 



u(z) dz 
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1 



2cos/?7rr(l - 2/3) 

a i 

2cos/37rr(l - 2/3) 

2£ 1 
2cos/3?rr (1-2/3) 

t 1 



— oo 
oo 



dz u(z) 

e l ^ z u(z) dz 
e^ z u(z) dz 



00 e**dx 



(x — z) 213 J-oo (z — a;) 



s l ^ x dx 



,26 



o y 



26 



dy 



o V 



26 



dy 



sin£y 



o V 



20 



dy 



cos /37r 


r(i- 


2/3) T(2/3) 




i 




cos /37T 


r(i- 


2/3) T(2/3) 


e 


i 




cos /37r 


r(i- 


2/3) T(2/3) 




i 


«(0 



/ / sin^ye-^u; 2 ^- 1 dwdy 
Jo Jo 



dw w 2 P 1 dye 



wy 



2i 



,2/3-1 



dw 



w 2 + £ 2 



cos/37rr(l - 2/3) T (2/3) 7 

_ e i g(0 r(/3)r(i-/3) _ 2/5 

2cos/37rr (1-2/3) T (2/3) |^| 2 - 2 ^ l?l ^' 

This concludes the proof of (3.13). 



(3.14) 



4. Space-time fractional telegraph equation 

We consider now the composition of an isotropic vector of stable processes S 2 @(t), 
t > 0, defined in (3.1), with the positively-valued process, defined in (2.19), 

L u (t) = inf is > : 9C{a) = H 2v {s) + (2\)^H%(s) >t\, t > 0, (4.1) 

where H 2v , H% are independent positively skewed stable processes of order 2v and 
i/, respectively. The distribution w@ (x,t) of the process S^f (c 2 X l/ (i)), t > 0, /3 g 
(0, 1], is the fundamental solution to the space-time fractional telegraph equation 



Cq2v 



dt 2 " 



m/3 



- + 2A— j w£ (x, t) = (- A) p ^ (*, t) 



x G 



,*>0. (4.2) 



a; e 



,t>0 



In our view the next theorem generalizes some previous results because we here 
have fractionality in space and time and the equation (4.2) is defined in R ra . 

Theorem 4.1. For v g (0, |], /3 G (0,1] and c > f/ie solution to the Cauchy 
problem for the space-time fractional n-dimensional telegraph equation 

((%£ + 2\&)w5(x,t) = -c 2 (-Af wt (x, t) , 

\wg(x,0) =6{x), 
coincides with the probability law of the vector process 

W n (t) = Sf (c 2 L»(t)) , t > 0, 
and has Fourier transform which reads 

S(€,f) = 



(4.3) 



(4.4) 
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_ 1 

~ 2 
where 



A 2 -c 2 ||£|| 2 ' 3 



E v ,i {nt") + l 



A 2 -c 2 ||£|| 2 ' 3 



(4.5) 



= -A + - c 2 ||£|| 2/3 , r 2 = — A — \/A 2 - c 2 |||" 2/3 



(4.6) 



T7ie time derivatives appearing in (4.3) mwst 6e meant in the Dzerbayshan-Caputo 
sense. The fractional Laplacian is defined in (3.6). 



Proof. By taking the Laplace transform of (4.3) we have 



H 2v wl (x,fi)- n^- 1 6(x)+2\ 



where we used the fact that (see [17] page 98, Lemma 2.24) 
-Cgu 

Now the Fourier transform of (4.7) yields 



C 



= -c 2 (-Afw v (x,n), 
(4.7) 

(4.8) 



( M 2 ^ + 2\(j, v ) wl (£, M ) - {l^- 1 + 2AM"" 1 ) = -c 2 \\tf 3 wl (€, . 
and thus 



At > 0,£ € R". 



(4.9) 



(4.10) 



( it 2!/ + 2A/z I ' + c 2 ||£|| 2 ' 3 ' 
The probability density of the process W n (t), t > 0, defined in (4.4), can be written 
as 

/>oo 

V£ (a;, c 2 s) C (s, i) <is, (4-11) 



(4.12) 



wP(x,t) = 

and has Fourier transform equal to 

J e* x < {x, t)dx = J e - c2s "«ll 2 " £,(«, t) ds. 

In order to show that the Laplace transform of (4.12) concides with (4.10), we have 
to derive the Laplace transform of 4,(x,i), with respect to the time t. Since 

Pr < x} = Pr{#"(x) > t} (4.13) 

we have that 
Z (x, M) = 



dx 



h v (s, x) ds I dt 



v— 1\ —xfj, u — 2Xxfi u 



(4.14) 



dx fi 

where we used result (2.28). Now we can complete the proof by taking the Laplace 
transform of (4.12) so that, in view of (4.14), we obtain 

/>oo />oo 

/ e'^dt / e-^ 3 ^ t v (s,t)ds = 
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o ~ m 2 - + 2A M - + c 2 ||€|| 2/3 ' 

(4.15) 

which coincides with (4.10). The unicity of Fourier-Laplace transform proves that 
the claimed result holds. The proof that the Fourier transform of (x, t) has the 
form (4.5) can be carried out by means of the calculation performed in Theorem 
2.2. We have that 



^ + 2A^ + c 2 ||||| 2 ' 9 M"-n ' n v — Ti {n v - n) (fi u - r 2 ) 



/x" — ri jti" — r 2 



2VA 2 -c a ||€|| v 

(4.16) 



where 



n = -A+ VA 2 - c 2 ||£|| 2/3 , r 2 = -A - a/A 2 - c? WtfP. (4.17) 



and thus by inverting (4.16) by means of (2.33), we obtain result (4.5). An alter- 
native derivation of (4.5) can be carried out as follows 

/>oo />oo 

«>£(€,*)= / e**-"dx Pr{S 2p (c 2 s) £dx}Pr{L u (t)eds} 

J-oo JO 

POO 

= e~ c2s ^ Pr{L v (t) € ds} = (4.5) (4.18) 
Jo 

because of Theorem 2.2. □ 

4.1. The case ^ = |, subordinator with drift. The fractional equation (4.2), 
for n = 1, v = h, reads 



\a + 2X m¥ ) w l {x > t] = c \WW ) ^ 4) ' ° < * < l! (4 - 19) 

where grL-gg is the Riesz operator defined in (3.12). For /3 = 1 we have the special 
case 

(l + 2A ¥r)^ (x ' t) = c2 £^ (x ' i) (4 - 20) 

dealt with in Orsingher and Beghin [21]. The construction of the composition 
related to equation (4.19) involves the subordinator 

= t + (2X) 2 H^(t), t>0, (4.21) 

where H 2 (i), t > 0, is a positively-skewed stable process and has the same law as 
the first-passage time of a Brownian motion through level We note that #"5 (t), 
t > 0, has distribution with support [t, oo) and thus differs from 9{ u (t), t > 0, 
< v < |, which instead has support [0, oo). The distribution of (4.21) writes 

x-t _j2 

Pr ksft <i = / -= ; rfz, x>i>0. (4.22) 
l > h \/2 v^FF 



22 



MIRKO D'OVIDIO, ENZO ORSINGHER & BRUNO TOALDO 



The inverse process 

£5 (t) = inf js : s + (2A) 2 (s) > t\ = inf |s : (s) > t\ (4.23) 
is related to (4.21) by means of the relationship 

Pr{xi(t) <af } = Pr {**(*)>*} = /" b (4-24) 

From (4.24) we can extract the distributon of £3 (t),t> 0, in the following manner 

Pr{x5(t) G dx} a ,00 £ 

4(M = " 3 ^ = "ET~ / ^ 

2 dx ax / t-* J/^^i 

(2A) a 

_ (2Aic) 2 _ (2Ax) 2 

2Aze *<*-<»> , e 4 <*— > 

2A 0<a;<t. (4.25) 



y^Tr ft - x) 3 



Remark 4.1. The distribution (4.25) can be also obtained from the general case 
(2.50) which for v = | becomes, for < x < t, 

k(x,t) = [ 8(s - x) hi (t - s, 2Xx) ds + 2A f h(2Xx, s) 8 (x - (t - s)) ds 
2 Jo 2 Jo 2 

= hi(t-x, 2\x) + 2A / 1 (2Ax, t - x) 

_ (2Ai) 2 _ (2Ax) 2 

2Aa;e 4 <*— 5 „, e 

2A - =. (4.26) 



In the last step we used the fact that 

Li(t) l = \B(t)\, t>0, (4.27) 

where (t), t > 0, dealt with in section 1.3, is the inverse of the totally positively- 
skewed stable process H2(t), t > 0. 

The t-Laplace transform of (4.25) becomes 

/>oo 

Tx(x,n) = / e" M * Ci (x, t) dt = 



X 



x/2 7 X ^27r(t-x) 3 7x - a:) 

V2 7o V2^rf3 7o V7rt 

= e -Mx e -2AxVM + 2A e -^ e -2AxVM^ (4 2 g) 

Finally the rr-Laplace transform of (4.28) becomes 
4 (7, A*) = jf e~ 7a; (J e'^ Ci(x,t)dt \ dx 

1 2A 1 1 + 2A//-2 



^ + 7 + 2A v ^ /i + 7 + 2Ay / /2 /x + 7 + 2A*//Z' 



(4.29) 
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which coincides with (2.31), for u = |. Let us now consider the process W n (t) = 
S?f (c 2 t > 0, dealt with in Theorem 4.1. For fj = 1, n = 1 and v = \ this 

process becomes 

Wi(t) = S 2 (c 2 X5(t)) = B (c 2 xs(t)) , i > (4.30) 

where B represents a standard Brownian motion and £a(t), t > 0, is the process 
defined in (4.23). With 

e 4 * 

P| B |(a;,i) = —j=i x>0,i>0, (4.31) 
V irt 



we denote the law of the process |-B(t)|, t > 0. In view of the previous results we 
are able to prove the following theorem. 

Theorem 4.2. The law of (4.30) coincides with the law of the composition 

W(t) = T (\B{t)\) , t>0, (4.32) 

where T is the telegraph process (1.26) with parameters c > 0, A > and law 
pr{x,t) which has characteristic function 



^ A | —\t-l-*. /\2_„2f2 I _ A 



-At+v^-^e 2 + i A ] e -xt-ty/\2- C 2e 

v/a 2 - c 2 e J " V va 2 - c 2 e ) 



(4.33) 

In other words we have the following equality in distribution 

B(c 2 L^{t)Y= T(\B(t)\), t>0. (4.34) 
Proof. First we show that the Fourier-Laplace transform of the law w\ (x, t) of the 

2 

process W\{t) = Sf (c 2 L^(t)j = B ^c 2 L^(t)j, t > 0, is written as in (4.15) for 
v = h, j3 = 1, n = 1, and reads 

g«,M)= ' + X 2 A^" 2 ,2 - ( 4 - 35 ) 



We have that 



3 />OG 

p B (a;, c 2 s) <is / e~^ L(s,t) dt 



poo poo 

/ p B (x,c 2 s)ds / i 



_ (2Aa)- i _ (2Aa)^ 

Alt , + 2A , dt 



OC 

PB 

n 

oo 



^47T (t - sf " vM* - S ) 

(x,c 2 s) (e-<»+ 2X ^) +2\JJle- s ^ +2X ^ ds 



/■mj Q — 

•('*+ a VP)da + 2A/i-a / ^L4 e - S (^+ 2A ^) ds 



e 4 C 2 S 

V 47TC 2 S J0 a/47TC 2 S 



(4.36) 
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and thus taking the Fourier transform we get 



w\ (£, fi) = 



Jo 



M + 2A^Z + c 2 C 2 ' (4 ' 37) 

Now we are going to prove that the law w(x, t) of the process W(t), t > 0, has 
Fourier-Laplace transform that coincides with (4.35). We have that 



(4.38) 



w{x,t) = I p T (x,s)p\ B \(s,t)ds, 
Jo 

and thus the Fourier transform of w(x,t) reads 
£>(&t) = 



1 

2 Jo 



e l ^ x dx I pt (x, s) p\B\{s,t) ds 
X 



o 



1 



VA 2 - c 2 £- 
A \ 



x v/a 2 - c 2 ," y 

Passing now to the Laplace transform we have 



e -As+s-vA 2 -c 2 £ 2 

P|B|(s,t) ds. 



As-s-^/A 2 -^ 2 



(4.39) 



M>(£,/i) 



2 Jo 



1 



A 



1 - 



A 







Jo 


[( 



v/A 2 - c 2 £ 2 

V 1 



e -A6- s% /A 2 -c 2 C 2 



-As+s-^/A 2 -^ 2 



e 4t 



7Tf 



VA 2 - c 2 £ 2 



-As+s\/A 2 -c 2 f 2 



1 - 



A 



2JJi 



^/a 2 - c 2 e 2 1 

A 



= -As- s% /A 2 -c 2 £ 2 



'VP 



1 



VA 2 - c 2 £ 2 / V A + \/M-\/A 2 -c 2 ^ 2 / 



1 - 



(a + v/a 2 - c 2 e 2 ) (a + + \A 2 - c 2 e 2 ) 

(2^V^ 2 - c 2 e) (M + 2A^m + c 2 £ 2 ) 

\A 2 - c 2 e 2 - a) (a + - \A 2 - c 2 c 2 



- c 2 C 2 j (/x + 2A^ A1 + c 2 £ 2 ) 
1 + 2A/x- 2 



H + 2XJJi + c 2 £ 



2<r2 ' 



(4.40) 
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which coincides with (4.35). 

This shows that for each t we have the following equality in distribution 



□ 



T(\B(t)\) l = B 



(c 2 X*(t)) 



t > 0, 



(4.41) 



where the role of the Brownian motion is interchanged in the two members of 
(4.41). Thus, by suitably slowing down the time in (4.41), we obtain the same 
distributional effect of a telegraph process taken at a Brownian time. 



Remark 4.2. The probability distribution of the process 



w\(x,t) — 



can be written as 

C7T 

A_ 

C7T 

y = \s VA 
C7T 



Wi(t) = B (c 2 L*(t) 



1 



t > 0, 



o ^/s(t- s) 
1 



e 4 ° 2y e 



e 4 = 



2(t - s) 

1 

2 



t 



t- s 

1 - 



ds 



ds 
t 

t - y 



Taking the limit for c — > oo, A — > oo, ^- — >• 1, formula (4.43) becomes 



(4.42) 



dy. (4.43) 



lim y\ (x, t) = 2 

A,c— >oo 2 



e 4 » e * 



/47TJ/ VTrt 



<7y 



(4.44) 



which coincides with the distribution of an iterated Brownian motion B\ {\B%(f)\), 
t > 0, with Bj,j = 1,2, independent Brownian motions. From (4.43) we can see 
that the distribution of W\{t), t > 0, has a bell-shaped structure. 

Finally we show that the density w\ (x, t) integrates to unity in force of the 

2 

calculation 



vu\ (x, t) dx 



dx 



e 4s 
ds ■ Ci (s. t) 
V4Vs 2 



ds 



d_ 
ds 



: dz 



se 4 * 

\/ 47TZ 3 



dz 



s=0 



te 4 * 

V 47TZ 3 



(4.45) 



In the intermediate step, formula (4.25) has been applied. 

Remark 4.3. The characteristic function of the process T 2 ^(t), t > 0, whose 
distribution satisfies 



'(j^ + 2Af) $(x,t) = c 2 £^p 2 T p (x,t), 0<p<\,P±\ 



p/(x,0) = 5 (a;), 



(4.46) 



at 



Pt (x, t) 



4=0 



0. 



reads 
Ee 
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1 + 



a/A 2 - c 2 |£| 2 £ 



o tV A2 - c2 i«i 2 ' 9 



+ 1- 



x /A 2 -c 2 iei 2 ' 3 



-tVA 2 -c 2 iei 2 ' 3 



(4.47) 



see Orsingher and Zhao [22]. Therefore by performing the same steps as in theorem (4.2) 
we prove that 

Sf (^(tj) '= T 2/3 (|B(t)|) , £ > 0. (4.48) 



4.2. The case f = ~, convolutions of Airy functions. We first recall that the 
totally positively-skewed stable process H^(t), t > has law 



Pr|if3(t) G dxj 



zAi 



x\/3x 



'3x 



dx, 



x > 0,t > 0, 



(4.49) 



where Ai(-) is the Airy function. Result (4.49) can be obtained from the general 
series expansion of the stable law of order | (see Orsingher and Beghin [24] page 
245) which reads 



k=0 



fc! 



3 ^-|(fc+D-l sin 



£(* + !>) 



. /27T(fe + l) 



3?r 

k=0 
1 3§ 



3 



Ai 



1 



'3x 



1 =Ai ( 1 



x V3x 



'3 a; 



where we used formula (4.10) of [24], which reads 

k sin 



Ai(tu) 



3 3 



>3 W 



27r(fc + l) 
3 



fe=0 



fc! 



fc + 1 



Since 



(4.50) 



(4.51) 



(4.52) 



we have result (4.49). From the relatioship between Hi(t), t > 0, and the inverse 
process Li (t), t > 0, 



Pr (t) < x j = Pr |i5 (a;) > i| 
we extract the density of Ls(x), x > 0, 



(4.53) 



Pr 



r |Ls (x) G rfij g ex ^ 



df 



1 



9t Jo s v^s 
1 = Ai 



zAi 



/3s 



ds 



o s V3s 



'3s 



ds- 



o sv 3s 



=Ai' 



t 



ds 



Since 



'3s/ v3s 

(4.54) 
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we conclude that, for x > 0, t > 0, 

Pr|is(s) e dt\ 
t 



ll (t, X 



dt 

-l 



'fix 



s\/3s 

-1 
s^3s 



=Ai 



Ai 



Ai 



/3s 



/3s 



ds 
(is 



3 d k .f t 
-Ai 



/3s 9s' 



/3s" 



,Ai 



/3s" 



3s. 

s—x 



ds 



r ds 
sMs" 



:Ai 



/3x 



i 

(4.56) 



In the last step we took into account the asymptotic expansion 7.2.19 of Bleistein 
and Handelsman [5]. 

With similar calculation we obtain the law hz{x,t) of the process H*(t), t > 0, 
which is expressed in terms of Airy function. From the general series expression of 
the stable law (see [24]) we have that, 

hg (x, 1) 



2 

3tt" 



2 ~ (-1)* g-f /fc + 1 



377^^ fc! 2 1 -K fe+1 ) 

v fc=0 



sm ( 



dwe u> 3 +2 



1 3 /22 1 r c 
x V 3a; 2 Jo 



e w w s Ai 



2 2 w 
3x2 



dio, 



(4.57) 



and thus, in force of the fact that H 3 (t) = t% Hz (1) 



hz(x,t) = 



t 1 r 



dwe w w a 3 



Ai -t 



22u 

3^ 



(4.58) 



Remark 4.4. We check that the distribution (4.58) integrates to unity. We have 
that 



hi (x, t) dx = 



dwe w w 6 \l — I dxx 3 1 Ai ( —t 



_ 1 3 npr 

y=x 3 t \ 



^ 3 



2 2 ( 



2 2 w 
3x2 
-1 



dy Ai(-y) 







1 

Jo 

1 

7^ 



1 3 / 2 2 / 3 / 2 2 



du> 1 



dwe w w 63 = 







dwe w w? 1 = 1, 



(4.59) 



28 



MIRKO D'OVIDIO, ENZO ORSINGHER & BRUNO TOALDO 



where we used the fact that 



dyAi(-y) = -. 

For the law of the process Li (x), x > 0, we therefore have that 
Pr jit (x) < t\ = Pr {#5 (t) > x\ 

ro ° ro ° t 1 3A22 



(4.60) 



Jx 



■ xz V : 2 Ai ( *V 3J 



e~ w w~^dwdz (4.61) 



and thus 



OO /'OO 



dWz 3/ 2 2 . / 3 2 2 w 



e 6 



OO /"OO 



Jx 



2 2 3 /2 2 



Z\/7r V 3z 2 V 32 



Ai' -t 



2 2 K 

Jz~ 2 



OO /"OO 



./a; 



Z\ it V 3z^ \ V i>z z 



e e 



OO pOO 



OO rOO 



Jo ' 3^ 2 

dwdz 3 [W ( 3 f¥w' 
V^V3^ I V 3^ 



1 9 4 . / 



22m 
3^ 



e •"w e 



00 f°° dwdz 3 [W A . / , ,/2*tiA 



2-y/7r Jo » 3a; 




e~ w w-e Ai -t 



2 2 u 
3^ 



dw. 



(4.62) 



For checking that (4.62) integrates to unity one can perform calculation similar to 
that of Remark 4.4. 

Now we have all the information to get the distribution of the process £ 3 (i), 
t > 0, by means of formula (2.50). We have that 



Ci(x,t) 



Pr 



dx 

lg(x,t — s)hi(s,2Ax)ds + 2A / li (2Xx, s) hg (t — s, x) ds 
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, - W -i 3 22 / J 2 2 w 

dwe w w e ? — — Ai 



V^(t -s) Jo y 3(t - s) 2 I y 3(t - s) 2 

2A f* , f°° , ,„ i./ J 2 2 u; \ . .f2Xx 



ds / dwe "°w eAi -i {7 — — A 



a: s 
2s + t - s 



3 s 2 2 
/3i V 3(* - s ) 



(4.63) 



Result (4.63) permits us to write explicitly the solution of the fractional telegraph 
equation (1.10) for v = |, f3 = 1 and n — 1, as 

2 

iUi(a;,t) = / Ci(s,t)ds, xeR,t>0. (4.64) 

4.3. The planar case. Let us consider the planar process 

T(t) = (X(t),Y(t)), t>0, (4.65) 

with infinite directions and finite velocity c, investigated in Orsingher and De Gre- 
gorio [23], which has probability law (see formula 1.2 therein) 

\ e -\t+±y/cH2-(x2+y2) 

r(x,y,t) = - ==-, x 2 + y 2 < c 2 t 2 ,t > 0, (4.66) 

K ' 2irc ^c 2 i 2 - (x 2 + y 2 ) V ' 

which satisfies the telegraph equation 

d 2 d\ ( d 2 d 2 \ 

- 2 + 2A- J r(x, y,t)=c 2 {^ + W2 ) r(x, y, t). (4.67) 

The distribution of T(i), t > 0, has a singular component uniformly distributed 
on the circle dC ct — {{x,y) £ K 2 : x 2 + y 2 = c 2 t 2 } with probability mass equal to 
e _At . The process T(t), t > 0, describes a random motion where directions change 
at Poisson paced times and the orientation of each segment of the sample paths is 
uniform in [0, 2ir). 

Let q(x,y,t) be the distribution obtained by means of the composition of the 
process T(t) with a reflecting Brownian motion with law 

_ s 2 

P\B\(s,t) = 6 ^L, <>0,s>0, (4.68) 
V irt 



which satisfies the equation 

-r4p|fl|(M) = (4.69) 
at 2 os 

and also 

d 2 

Q- t P\B\(^t) = ^P|fl|(M) ( 4 - 7 °) 
We have the following theorem. 
Theorem 4.3. The law of the composition 

Q(t) = T(\B(t)\), t>0 (4.71) 
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written as 



Q{x,y,t) 



r(x,y,s)p\ B \(s,t) ds, 



satisfies the 2- dimensional time-fractional equation 



d c <93 

dt + 2X ^t¥ 



q(x,y,t) = c 

subject to the initial condition 



\ dx 2 dy 2 



q(x,y,t) , 



q(x,y,0) = S(x,y). 
Proof. By considering (4.72) and (4.69) we can write 



(4.72) 

x,y e R,t > 0, 

(4.73) 

(4.74) 



-z-r<l(x,y,t) 
at 2 



r{x,y,s) — —p\B\(s,t) ds 
at 2 



( d \ 
r{x,y,s) \ -—p\ B \{s,t) J ds 

P\B\{s,t)r(x,y,s)] s s J^ + p\ B \(s,t)—r(x,y,s)ds. (4.75) 



In the previous step it must be taken into account that the boundary dC cs is 
excluded. From (4.72) and (4.70) we have that 



d 

at 



q(x,y,t) 



d 



r(x, y, s) -Q t P\B\ (s, t) ds = J r(x, y, s) q^P\b\ (s, t) ds 

d .a . . , 

g^nx,y,s) -g s P\B\(s,t) ds 



a 



r{x,y,s) g^p\ B \(s,t) 



P\B\(s,t) ^r(x,y,s) 



•'O 

s—oo 



=0 Jo 



d 2 

P\B\{s,t) -^r(x,y,s) ds. 



(4.76) 



Thus, by looking at (4.67), (4.75) and (4.76) we obtain 

d c d^ 
w,q(x,y,t) + 2X— 1 -q(x,y,t) = 

at dt 2 



d 2 d 
Q-2 r (x, y, s) + 2\ — r(x, y, s) 



P\B\(s,t) 

f°° ( d 2 d 2 \ 

' P\B\(^)^[g^ + Q^)r(x,y, S )d 



ds 



\ dx 2 dy 2 



which means that q(x,y,t) satisfies equation (4.73). 



q(x,y,t). 

(4.77) 
□ 



It is easy to show that the process Q(t) = T (\B(t)\), t > 0, has not the same law 
of the process W2(t) = B2 {^c 2 L^(t)^j, t > 0. However it is possible to construct a 
planar process, say X(f), t > (which is a slightly different version of T(t), t > 0) 
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composed with a suitable "time process" which has the same distribution as W2(t), 
t > 0. The planar random motion %{t), t > 0, with distribution 



t(x,y,t) 



Ae 



-A/ 



2vrc 



; A^ c 2 t 2_ (a ,2 +y 2 ) + e _A^ c 2 t 2_( a ,2 +2/ 2 ) 

v/c 2 * 2 - (x 2 +Fy 



(4.78) 



where (a;, y) G C c t = {(a;, y) : x 2 + y 2 < c 2 t 2 }, can be constructed starting from the 
model dealt with in Orsingher and De Gregorio [23]. The distribution is based on 
the solution to the planar telegraph equation 

, ( d 2 d 2 

t(x,y,t) = c z 



d 2 d 

W 2+2X dt 



dx 2 dy 2 



*{x,y,t), 



namely 
v(x,y,t) = 



,-At 



^/c 2 t 2 - (x 2 + y 2 ) 



with A = B = and thus we can easily check that 



dxdyx{x,y,t) = 1 — e 



-2At 



(4.79) 



(4.80) 



(4.81) 



c ct 



We take a particle starting from the origin, moving at finite velocity c, and changing 
direction (chosen with uniform distribution) at Poisson times and neglect displace- 
ments started off by even-labelled times. The sample paths of this motion are 
constructed by piecing together only odd-order displacements of the planar motion 
T(t), t > 0. The process just described has distribution (4.78) as shown below 



t(x,y,t) = 

Pr {%(t) G dx} 



X 2 1 



dx 



\e- xt 
2ttc 



.fc=o 



J2(-^/c 2 t 2 -(x 2 +y 2 ] 



y/c 2 t 2 - (x 2 + y 2 ) 

2fe-l , 



^E(^) (2k + l)(c 2 t 2 -(x 2 + y 2 )) 
k=o v J 

oc 

2 Pr {X(t) G da?, F(t) G dy|iV(t) = 2fe + 1} 



(2fc)! 



,-Ai 



(\t) 2k 



2k+l 



k=0 

oo 



(2fc)!(2fc + 1) (Xt) 2 ^ 1 

-At (At) 2fc+1 
(2fc + l)! 



= 2 ^ Pr {T(i) G dsc|iV(t) = 2fc + 1} 



2fc+l 



fe=0 



(2fc + l)!' 



(4.82) 



where, for x 2 + y 2 < c 2 t 2 (see [23]), 
Pr{X(i) G dx,Y(t) G dj/|iV(t) = n} 



and 



fc=0 



dx dy 

{\t) 2k+1 
(2fc + l)! 



2n(ci) 

OO 

^ 2 Pr{ AT(t) = 2fc + l} = 1-e 



(^-(x^y 2 ))^ 1 , (4.83) 



•2At 



(4.84) 
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The factor 2 appearing in (4.82) and (4.84) can be interpreted as follows. The 
displacements generated by an even number of Poisson events are disregarded and 
replaced by displacements produced by an odd number of deviations. Therefore, 
odd-order Poisson events ignite twice the displacements considered in (4.82). 



Theorem 4.4. The composition with distribution 



<\(x,y,t) = / dsv(x,y,s) 



1 <95 

P\B\ (S,t) + 7y^-^TP\B\ (S,t) 



\ dx 2 dy 



+ 2 \— r } q (x,y,t) = c 2 [— + — )^x,y,t), 



which satisfies the time-fractional equation 
( d °d*' 

has the same law of the process W^i) = -B2 (^c 2 L^(t) 
Proof. We begin by evaluating the Fourier-Laplace transform of (4.85). 



(4.85) 



(4.86) 



OO fOO 



JO 

2\ + ^Jl 



ds / die-"* / dxdye i(x+iay t(x,y,s) 



( '< • 



1 95 

P\B\{s,t) + — — jp ]B] (s,t) 



2A^/m Jo 



ds I dx dy e iix+iay r (x, y, s) e~ s ^. 

Get 



(4.87) 



Now we need the Fourier transform of the law t(x, y,t) of the process X(t), t > 0, 
which reads 



t(£, a, t) 

Xer xt 
2-kc 



i£x-\-iay 



Cat 



yJcH 2 - (a; 2 + y 2 ) 



dx dy 



\ -At r 27T r ct \ P ^V c2t2 -P 2 -l P -T , A.2+2 _ n 2 

^— / de / d PP e^ cose+asi ^-— +e Vct P - 

2ttc Jo Jo c 

2 \2 -At fCt °o / ^ s, 2m-l 



JU rn=0 V 



(2m) 



^/c 2 * 2 - p 2 

^ Jo (We 2 + « 2 ) ^ 



2Ae^ ~ (J£ - 
c ^ (2m)! ^ 

m=0 v ' k=0 



(-1)* (^) 



\Y2 



2 k 



.... v ^ / v k a 



(c 2 t 2 - p 2 )" ^ P 2k+1 dp 



, 00 / , \ 2m 

m=0 v 



' Q (2m)! ^ 2(fc!) 2 (c*)-( 2 "»+ 2fc + 1 ) J 

^y-to } \ 2 ) 



y«{l-y) m -*dy 



2k 



{ct) 2m+2k+1 T(m+\) 



k\T (k + m + 1 + 



(4.88) 



SPACE-TIME FRACTIONAL TELEGRAPH EQUATIONS 



33 



Thus, from (4.87), we have that 
2X + y /Ji 



^' a ' M) 2AVM Jo 
in force of the calculation 



1 + 2X/j,- 



fi + 2X^I + c 2 (e + a 2 ) 



(4.89) 



POO 

/ dst(£,a,s) e~ 
Jo 



„oo oo , 2m oo 

1 dSe ' XS ^c^{2 m )\^ 



{-If 



2 



2/v 



r(m + |) 



^ o c 2 ™(2m)! ^ fc! ( cs )-(2m+2fc+i) r (fc + m + 1 + |) 



(-1)* ( ^ 



2 A 



„2A' 



_ y, 0f2 1 - 2m r(2m) 

ho A- 2m (2rn)!r(m) ^ fc!r (fc + m + 1 + §) 7 



E 



2 



2k 



T (2fc + 2m + 2) 



2 (A + Jjif to ml (A + ^) 2m ^ fc! (A + ^i) 2k r (fc + m + 1 + \) 



"A g A 2„- 2 l-2„, 



E 



(-i) A 



2A- 



r (fc + m + 1) 



2 (A + VM) 2 i=o ml (A + m to k\ (A + ^) 9 * c"" 2i- 2 (*+™+D^ 



2A 



E 



E 



. 2fc 



,-2/c 



2A 



(A + ^) 7o 



A 2 c 2 Q 2 +<» 2 ) 

u 7 rTj- — U — j-^ 

due 0+^) ( A +^) 



2A 

(a+7m) 2 



, c^(C 2 +a 2 ) 



(A+yp) (A+7m) 



2A 



2A 



(A + y/p) 2 - A 2 + c 2 (£ 2 + a 2 ) M + 2A^ + c 2 (£ 2 + a 2 ) ' 



(4.90) 



The Fourier-Laplace transform of the law of the process [c 2 L^ (t)j is written 
as in (4.15) for n = 2, P = 1 and z/ = | as the following calculation shows 



w\(£,a,t) = / vb (£,a,c 2 s) Ci (s,n) ds 



1 + 2A/i 



e -2As^ + 2 A 



2XJjl + n + c 2 (e+a 2 ) 



ds 

(4.91) 



In the previous calculation we use the Laplace transform of Ci(x,t) obtained in 
(4.28). The proof is complete since (4.91), coincides with (4.89) and with the 
Fourier-Laplace transform of (4.86). □ 
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Remark 4.5. Since for the first passage time = inf |z : B(z) = ^75} of a 
Brownian motion through level ^ we have that 

J e^'Pr |r_^ e dt} = e" s ^, (4.92) 

/ e-"*-rP| B | = e- s ^ (4.93) 

Jo ot 2 



and 



we can write 



t(x,y,s)—TrP\ B \{s,t) ds 

2 

d i \ e ~^ a 
^-*{x,y,s)—=ds = 

OS J-Kt 



d_ 
ds 



t(x,y,a)p\ B \(a,t)ds. (4.94) 



This representation of the second term of (4.85) is extremely interesting because 
by integrating (4.94) in C ct we get 

d 



(1 - e- 2Xs )p lB{ (s, t) ds = 2A / e-' As p lB] (s, t) ds (4.95) 



-2As 



/o ds 

and yields the missing probability of the first term of (4.85). 
Remark 4.6. We check that the law 



1 95 

P\B\{s,t) + ^x^tP\b\( s ^) 



ds (4.96) 



integrates to unity. By taking the i-Laplace transform, the integral with respect to 
(x,y) becomes 

/>oo 

dxdy / dte~ lit q(x,y,t) 



C c 

OO 



1-e- 



1-e 



-2As> 



-2 As 



e p ( I P| B |(s,i) + — ^rP| B |(s,t) I 



c/.s 



1 | 1 
+ 2A 
2A + / 1 



e- s ^ds 
1 



e -,(2A +y7 z) ds 

/•OO 

e'^dt. 



, _ . _ _,=-=/ e-"*dt. (4.97) 

The same check can be done directly by taking into account formulas (4.94) and 
(4.95). 

Relationships similar to B (c 2 L^ (t)\ ^= T (\B(t)\), t > 0, and the analogous one 

in the plane, cannot be established in spaces of dimension n > 3, because random 
motions governed by telegraph equations in such spaces have not been constructed. 
Random flights in M. n have been studied (Orsingher and De Gregorio [23]) but their 
distributions are not related to higher-dimensional telegraph equations. 



SPACE-TIME FRACTIONAL TELEGRAPH EQUATIONS 



35 



5. Hyperbolic fractional telegraph equations 
The Hyperbolic Brownian motion is a diffusion on the Poincare half-space 

H" = {(x,y) :xGM n -\y>0}, (5.1) 
with generator, written in cartesian coordinates, 



fin = 



1 



d 2 



d_ 

dy 



(5.2) 



In the half-plane H 2 the hyperbolic Brownian motion was introduced by Gertsen- 
shtein and Vasiliev [12] while in H 3 it was introduced by Karpelevich, Tutubalin 
and Shur [16]. The reader can also consult, for more details, Getoor [13]; Gruet 
[14]; Lao and Orsingher [18]; Matsumoto and Yor [19]. The hyperbolic Poincare 
half-space is equipped with the metric 



ds 2 



y 2 



(5.3) 



and thus the hyperbolic distance in HP is given by the formula 

cosh^z'.z) = 1+ ^ Z '~ Z } , z,z'eH n , (5.4) 

where ||-|| is the usual euclidean norm. We define the operator fi 2 &s the governing 
operator of the planar hyperbolic Brownian motion B2 P (t), t > 0, which is written 

as 



S)2 



\ dx 2 dy 2 



in Cartesian coordinates and takes the form 



in hyperbolic coordinates, where 



Q2 



1 



d 2 



sinh 2 r\ da 2 



d 



sinh i] drj 



sinh?7 



(5.5) 



(5.6) 



(5.7) 



Note that we disregard the factor | in Sj 2 P m the forthcoming calculation as in 
the pioneering work by Gertsenshtein and Vasiliev [12]. The problem involving the 
radial part of (5.6) which is written as 



7] > 0,<> 0, 



\k 2 (ri,0) = %), 
has the following solution 

_ t 

k 2 (ri,t) 



23 V / 7ri 3 



ipe 4 * 
y/ cosh ip — cosh r\ 
hp i 



dtp 



(5.8) 



(5.9) 



to which we refer as the kernel of the law of B% p (t), t > 0. The law of B% p (t), t > 
is therefore written as 

P2 P (Vit) — sinhry ^2(7/, t), 77 > 0,< > 0. 



(5.10) 
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The three-dimensional hyperbolic Brownian motion B^ p (t), t > is driven by the 
operator 

i0 3 = z 2 (— + — + — )-z— (5 11) 

\dx 2 dy 2 dz 2 J dz 

written in Cartesian coordinates. We are interested in the Cauchy problem 
fw k z(V,t) = gsfta (»?,*), r/>0,t>0, 



Afcfa.O) = 5(7?). 



where 



represents the radial part of Sj^ p which coincides with Sj 3 in hyperbolic coordinates. 
The solution to (5.12) is given by 

-t -it 

k 3 (r,,t) = e **\ (5.14) 
2\/7rf3 sinhTj' 

and thus the probability law of B^ p (t), t > 0, reads 

p% p (r),t) = S mh 2 r}k 3 (r),t). (5.15) 

In general, the law of a n-dimensional hyperbolic Brownian motion is written as 

d*(r,,t) = smbT-^hnfat), (5.16) 

and solves the heat equation 

where 

a-sK-^'IKsiM)- nEN ' (5 ' 18) 



is the adjoint of 



in the sense that 



(g n k n , Pn) = {K, GnPn) , n G N. (5.20) 

Thus the n-dimensional kernel satisfies 

—k n (r),t) = g n k n (f],t). (5.21) 

The kernels for n > 3 can be obtained from k% and fea by means of Millson 
recursive formula (see Debiard, Gaveau and Mazet [7]) 

e -ni Q 

fcn+2 (*?,*) = ~ ■ r (*?.*)■ (5-22) 
2tt smh ry #77 
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By working out the derivatives we obtain a more explicit version of Millson formula 



~V 1 ) t ( a V 



> 1, n = 2j + 1, 



k 2j+2 {r),t) = e {27[}] 



h (r],t) , 



Theorem 5.1. The distribution of the composition 
<tf(t) = B%>(L"(t)), ve(o 



j>0,n = 2j + 2. 



t > 0, 



(5.23) 



where B^ p is the n- dimensional hyperbolic Brownian motion in the Poincare hy- 
perbolic half-space H n , satisfies the fractional hyperbolic telegraph equation for v g 

UM) = S(r,), 
and thus the kernel 



satisfies, for v G (0, |] . 



sinn 77 



(5.24) 



( 



"»"'+ 2A- 



at2 „ -r *n gt 



(5.25) 



Proof. It is convenient to consider the Laplace transform, of the kernel K v n . We 
have that 



dte-^K^fat) = / dte-^ dsk n (r},s) ^(s,t) (5.26) 
Jo Jo 

k n ( v , s) (fi 2 "- 1 + 2X^- 1 ) e " V+ 2A ^) ds. (5.27) 



Now we show that (5.27) satisfies the Laplace transform of (5.25) written as 

1 



(^ + 2A^)< (r?, M ) 



sinh" 1 rjdr) 



d ( . __i d ~ . 



(5.28) 



By considering (5.21) and that (77, 0) = 5 (77) we have, for 77 > 



1 g 
sinh" -1 77 <9?7 



fsinh" 1 77^^(77, /j ) 



77 y 

* * (sinh- 1 ^M*-)) (^- 1 + 2A M - 1 ) e -(p*W) 



sinh" 77 

ds^-k n ( V ,s) ( / x 2 "- 1 +2A/i"- 1 ) e -(^W) 
as 



s— 00 
s=0 



= fc„(77, s) (/j 2 - 1 + 2A / i- 1 ) e- s ("" +2A "") 

/•OO 

+ (/j 2 ^ 1 + 2A^- 1 ) (/a 2l/ + 2A/i") / fc n (77, s) e - ds 

Jo 



(5.29) 
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□ 



Remark 5.1. By taking profit of the simple structure of p^ p (77, t) we can give, for 
n = 3, an alternative direct proof of the result of theorem 5.1. We first evaluate 
the Laplace transform K^(rj,fj,), as follows 



K 3 (V,!*) 



fc 3 (?7,s) / e ^(visttydtds = / k 3 (rj,s) C u (s,[±) ds 



-s(l+^ 2 "+2A^) 



2y/TT sinh 77 



ds 



sinh 7] 



(5.30) 



Now we show that (5.30) solves the Laplace transform of (5.25) for n = 3. We have 
that 



1 d 



I 2 d ~ 
( sinh J7^«3 (r?,/i) 



a (z^- 1 + 2A„»-i) , n/ l+ ^ +2AM , 



sinh 2 77 9?7 V 

sinh 2 77 9ry I <9?7 sinh 77 

_ (/J +|V ) 9 L-r,Vl+M 2 "+2A Al - [ sinhr/ ^1 + fj2u + 2Xfj," + cosh 7/) 

sinh rj or] L V / 

T-T2 V 1 

sinn 77 l 
• ( sinh 77-^ 1 + il 2v 



+ ]i 2v + 2Xfi"- 
2\\i v + cosh 77 j - (cosh 77^1 + \i 2v + 2\\i v + sinh 77 

^ P : % [sinh 77 (1 + /i 2l/ + 2A/i ,y ) - sinh 77] 

sinh 77 

(/i 2 "" 1 + 2A/J- 1 



sinh?7 



(/i 2l/ + 2A/i 



(5.31) 



Remark 5.2. For ^ = g we know the explicit law of the process L v (t), t > 0, 
which is written as in (4.25). Thus we have an explicit representation for the law 
of the process 



(5.32) 



which reads 

Pi(v,t) = sinh 2 77 



?7 e 



2\/7rs3 sinh 77 
A 77 sinh?7 



Ase 



2Ae~ 



2tt 



S2 ^/t — i 



\/7r(i™s)3 \Ar(i - s) 
— ^— + 2 1 ds. 

E — S 



r/.s 



(5.33) 
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